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Abstract 

Rho meson propagation in hot hadronic matter is studied in a model with 

coupling to vrvr states. Medium modifications are induced by a change of 

the pion dispersion relation through collisions with nucleons and A's in the 

fireball. Maintaining gauge invariance dilepton production is calculated and 

compared to the recent data of the CERES collaboration in central S+Au 

collisions at 200 GeV/u. The observed enhancement of the rate below the rho 

meson mass can be largely accounted for. 
21.65.+f, 25.70.Np. 
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Typeset using REVT^ 



In central heavy-ion reactions at ultrarelativistic energies hot and dense hadronic matter 
is formed at the early stages of the collision. From such experiments one hopes to infer 
properties of mesons and baryons in the vicinity of the chiral phase transition. Because 
of its electromagnetic origin, lepton pair emission is an ideal probe. Once produced, the 
dilepton will decouple from the strongly interacting particles and thus carry undistorted 
information of the dynamical properties of the hot nuclear system. For this reason the 
recent low-mass dilepton spectra measured by the CERES collaboration in p-|-Be, p-l-Au and 
S-|-Au IQ collisions at 450 GeV/u and 200 GeV/u, respectively, have produced considerable 
interest. While the proton data are well described by known decays from a hadronic fireball 
an excess of dielectrons with invariant masses of 200 MeV to 1 GeV has been observed in 
the S-|-Au case and has been attributed to tt+tt^ annihilation from the pionic component of 
the hadron gas . When this process is included in the description of the data the excess of 
strength between the 2m7r-vacuum threshold and the rho peak is still unexplained. This led 
the authors of ref. to the conclusion that one is seeing a signal of partial restoration of 
chiral symmetry through a decrease of vector meson masses (here especially the rho meson) 
near the phase boundary. Since the rho meson strongly couples to the tt+tt^ channel this 
effect may also be explained more conventionally, however, namely by modifications of the 
pion propagation in the presence of nucleons and isobars f^-^. 

To study the dilepton emission from tt~^7c~ annihilation we invoke the vector-dominance 
model (VDM) in which the electromagnetic current is related to the third isospin component 
of the rho meson field by the field current identity J'^ = {m?p/ g)p'^ where and g denote 
the rho meson mass and the universal VDM coupling constant, respectively. In terms of the 
retarded (thermal) rho-meson propagator 

r p-f3Ei 




(1) 



(/3 = 1/T) the dilepton emission rate per unit volume, R=dNi+i~ / d^x, is given by 



H{q-T) 



(2) 
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where q = p+ + P- = {(lo,(l) denotes the total momentum for (massless) dileptons, = 
> is the invariant mass and the hadronic piece H{q] T) is given by 

H{q;T) = nqo,T){imf)y7rg')g,amD^/{q;T) (3) 

(f^(qQ,T) = (e'"'^'^ — 1)^^). To reduce the numerical complexity in evaluating D'^" we shall 
restrict ourselves to back-to-back kinematics, i.e. q = 0. In that case, from gauge invariance, 
only the space components contribute. Denoting D^J{qo, q = 0;T) = 6^^Dp{qo; T) eq. (3) then 
reduces to 

H{q,,q = 0;T) = -/^(go, r)^^ImD,(go; T) (4) 

where 

DMo: T) = [ql - {mff - S,(go; T)]-' (5) 

involving the scalar part of the rho self energy Sp. 

To study medium modifications of the p propagator we first need a realistic model for the 



p meson in free space. As is well known |T0| an appropriate description is provided by a bare 
pole graph renormalized by tttt rescattering. This amounts to solving a Lippmann-Schwinger 
type equation 

oo 

M,,(E, gi, gs) = Kr.(^, gi, 92) + J dk ^ V^^{E, gi, k) Gl^{E, k) M^^{E, k, q^) , (6) 



where E is the starting energy of the pion pair, gi and g2 denote the moduli of the pion 
three-momenta in the cms frame and G^^{E, k) is the free two-pion propagator: 

GlAE,k) = - } ■ ul = ml + k'. (7) 

ujk E^ - + IT] 

The pseudopotential generated by the s-channel p pole graph is given by 

V^^{E,q^,q2)=v{q^)Dl{E)v{q2) , (8) 

where 



are the bare p propagator and the vertex functions, respectively. The hadronic form factor is 
chosen to be of dipole form normahzed to one at the physical resonance energy mp=770 MeV: 

- omt ■ 

The parameters g, Ap and m^°^ are fitted to the experimental p-wave tttt phase shift and the 
pion electromagnetic form factor 

(^2 - (mf )2 - /?eS0(E))2 + JmS0(E)2 
^ (mf^r \DliE)\' (11) 

with 

s°(i?) = / edkv{kyGi{E,k). (12) 

The subtraction for S° at E=0 is necessary to ensure the normalization -^,^(0) = 1. We 
obtain a satisfactory fit with g"^ /A7[=2.7, Ap=3.1 GeV and m^p^=829 MeV as indicated in 
fig. 1. 

A proper calculation of the in-medium p propagator 

Dp{E) = [E'-{m'-^^r-i:p{Er\ (13) 

must include constraints from gauge invariance. A detailed derivation for the self energy in 
homogeneous matter is given in refs. |^^, and we do not repeat it here. The result can be 
expressed in terms of the spin-isospin longitudinal and transverse response functions at a 
given density and temperature 

U^iko, k) = {kl -e- ml) n°(A;o, k) DL{ko, k) , 

Uriko, k) = {kl -e- ml) ti^k^, k) Driko, k) (14) 
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where Dl and Dt are the longitudinal and transverse propagators which are expressed in 
terms of the pion self energy as 

Di(fco, k) = [kl -e-ml- S.(A;o, fc)]"' , 



Dt(A:o, k) = [kl -e-ml- Cp^Ako, k)]-' 
1 
k 



ii\ko,k) = -j:^{ko,k) . (15) 



Obviously -Dl is equal to single pion propagator D^^. Making use of derivations in ref. 
the imaginary part of the self energy for a p meson at rest then takes the form 



ImSp(go, 0) = J k^ dk vik)^ J — ^ImD^(go — ko, k) 





xIm{a(go, ko, k)D^{ko, k) + ]^TiL{ko, k) + IiT{ko, k)], (16) 

where the function a is given by 

a(go, A;o, A;) = [1 + n°j(A;o, k) + Iil{qo - ko, k) + Uil{ko, k)Iil{qo - ko, k)] 
fl%{ko,k) = Refl''{ko,k) . (17) 

The real part of is obtained via a dispersion integral: 



The factor Qq/E'"^ results from a subtraction at qo=0, which is necessary to ensure gauge 
invariance f^; the imaginary part is unaffected by this subtraction since /mSp(go = 0) = 0. 

It remains to specify the model for the single-pion self energy S^. We employ the 
standard model of p-wave particle-hole excitations extended to finite temperature P]. Here 
the pions are dressed by NN^'^, AN^'^ as well as A^A"-*^ and AA^^ excitations, the latter 
appearing as a consequence of a thermally excited A abundance in the gas. One has 

J:^{uj,k) = -e J2Xa{oo,k) , (19) 

a 

where the summation is performed over all excitation channels a=ab~^, a,b=N,A. The 
susceptibilities Xa contain short-range correlations taken into account by Landau-Migdal 
parameters g'^/^ such that 
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Y _ ^(0) _ (0) / 

X^^\uj,k) = (IjUllMy SI {a) U^^k) , (20) 
\ ) 

where f-j^a denote the irNN etc couphng constants and SI{a) is a spin-isospin factor. The 
form factor T.,^{k) = (A^ — m^)/(A^ + k"^) accounts for the hadronic size of the pion-baryon 
vertex (A^r = 1200 MeV). The exphcit form of the thermal Lindhard functions, (pa, including 
the A width has been given in ref. 0. For the Migdal parameters we take g'^^ = 0.8 for 
a(5=aa~^hh~^ and g' = 0.5 for all others. Finally, to account for a finite pion density, 
we supplement eq. (16) with a 2-pion Bose factor, [1 + f^iko) + /'"(go — ^o)]? which can be 
derived rigorously within a Matsubara formalism [|TT|. The contribution of the gas pions to 



the pion self energy ^^^(/co, k) has been shown to be small [TT| and we can safely neglect it 
here. 

Expressions (2)-(20) specify our model for the dilepton rate in back-to back kinematics. 
Since we evaluate the p propagator for a homogeneous gas of nucleons and A's in thermal 
equilibrium, the rate at a given temperature T of the fireball is determined by the nucleon 
and A abundances at that T. To evaluate these abundances we make use of the transport 
results of Li et al. p|. It is found that at the initial stage the total baryon density pb is 
3.5po (Po = 0.16 fm~^) at a temperature of 170 MeV. Including besides the nucleon and A 
all baryon resonances with masses below 1.7 GeV as well as the lowest-lying hyperons 
and assuming chemical equilibrium, the chemical potential p^ = p/^ = . . . = 0.448 GeV is 
determined by the initial baryon density. The corresponding nucleon and A densities are 
1.0 Po each, the sum being in good agreement with the results quoted in ref. 0. The time 
evolution of the abundances can be obtained from the time dependence of the temperature 
given in ref. 0. A reasonable fit of the transport results is obtained with 

T{t) = {T - T°°)e-*/^ + T°° (21) 

with an initial temperature T*=170 MeV, T°°=115 MeV and r=10 fm/c. For chemical equi- 
librium the time evolution of pb can then be calculated and the result is again in agreement 
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with the transport calculations p[. Thus we can extract pn and pA at each time for the 
integration of the time history. 

For a direct comparison with experiment two further points have to be considered. The 
first is related to the fact that the rate is only evaluated in back-to-back kinematics (^=0). 
Minimally this can be corrected for by using []T2| 

dNin- ( dNi+i- \ 
F(M,,,T)^-P[-(-f^-^)/^] (22) 



+ 

which is exact for free pions. The second, more severe, point is the finite momentum accep- 
tance of the detector. In the CERES S+Au experiment only dielectrons with opening angles 
6ee > 35 mrad and transverse momenta pt > 0.2GeV are detected. While the opening angle 
restriction is not severe the pt-cuts have to be included properly. To do so we statistically 
model e~^e~ p decays at given temperature and invariant mass and apply the acceptance 



cuts to each of the two lepton tracks |jT3| resulting in an acceptance function A{M,T). The 
final dilepton yield can now be calculated by taking into account the above modifications 
and integrating over the time evolution of the fireball up to the freeze-out time to obtain 



d^e+e- _ J^, ^f^, rT.u^^ fdNe+e-it)\ 



where 



C(M, T) = A{M, T) J ^F{M, g, T) (24) 

contains the acceptance function as well as the finite q correction to the back-to-back rate 
(eq. (2)). We take t/=10 fm/c, but our results are affected less than 10% when choosing 
t/=15fm/c. Finally we fix the overall normalization in the measured rapidity window to 
the total yield from vacuum p decay obtained by Li et al. |^ . Fig. 2 shows our final results 
supplemented with contributions from free Dalitz decay of vr^'s, r^'s and cu's (dashed-dotted 
line), which we extracted from ref. We also implemented the contributions from free uj 



7 



decay {u e+e~) according to ref. 0, which, due to the small uj width, are expected to 
undergo only minor medium modifications 0. The inclusion of the medium modifications 
of the p meson leads to an enhancement of the e+e^ yield of about a factor of 3 for invariant 
masses around GeV/c^ (full line compared to the dotted line in fig. 2). 

In summary, we have presented a gauge invariant calculation of the p propagator at rest in 
a hot ttNA gas taking into account the full off-shell dynamics of the intermediate pions when 
coupled to nucleons and A's. After correcting for a finite momentum of the dilepton pair 
as well as for the experimental acceptance of the CERES detector, our calculated yield can, 
to a large extent, account for the enhancement in the e^e^ spectra observed in 200 GeV/u 
S+Au collisions in the CERES experiment. 
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Figure Captions 



Figure 1: Our fit to tlie p-wave vrvr pliase sliifts (upper panel) and tlie pion electromagnetic form 
factor (lower panel); the squares in the lower panel are the values from the Gounaris- 
Sakurai formula [|14[, which itself gives an accurate description of the data. 

Figure 2: Dielectron yield from free Dalitz decay (dashed-dotted line), free Dalitz + free uj + 
free p (dotted line) and free Dalitz + free ui + in-medium p decay (full hne); the dots 
are the CERES data. 
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